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Abstract. We study the effect of spin-phonon coupling in small magnetic clusters, concentrating on a 
S=l/2 ring of 4 spins coupled antiferromagnetically. If the phonons are treated as classical variables, there 
is a critical value of the spin-phonon coupling above which a static distortion occurs. This is a good 
approximation if the zero point energy is small compared to the energy gain due to the distortion, which 
is true for large exchange interactions compared to the phonons energy (J > huj). In the opposite limit, 
one can integrate out the phonon degrees of freedom and get an effective spin hamiltonian. Using exact 
diagonalizations to include the quantum nature of both spins and phonons, we obtain the spectrum in the 
whole range of parameters and explicit the crossover between the classical and quantum regimes. We then 
establish quantitatively the limits of validity of two widely used approaches (one in the quantum and one 
in the classical limits) and show that they are quite poor for small magnetic clusters. We also show that 
upon reducing Hu)/J the first excitation of a 4-site cluster becomes a singlet, a result that could be relevant 
for CuaTeaOsBrz. 

PACS. 70 Condensed matter: electronic structure, electrical, magnetic, and optical properties - 75 Mag- 
netic properties and materials - 75.75.+a Magnetic properties of nanostructures - 75.45.+j Macroscopic 
quantum phenomena in magnetic systems 



1 introduction 

Although electron-phonon and spin-phonon couplings are 
always present in all materials, at low temperature the 
phonons are usually frozen and many properties of solids 
are correctly described by pure electrons and/or spins 
hamiltonians, as for instance the Heisenberg hamiltonian. 
There are however some examples where this coupling 
■ is primarily responsible for the low temperature physical 
properties. Superconductivity is the most famous example 
among electronic systems. In magnetic systems, the spin- 
Peierls instability in the antiferromagnetic spin | chain 
and also in some frustrated two dimensionnal systems 
[2] shows the importance that spin-phonon coupling might 
play. These instabilities with respect to a spontaneous dis- 
tortion of the lattice for any value of the spin-phonon cou- 
pling are special cases, and in standard systems there is 
a finite critical value of the spin-phonon coupling above 
which the system will be distorted. It is in particular the 
case for finite systems such as rings consisting of a small 
number of spins [3]. In the case of a degenerate ground 
state, a distortion lifting this degeneracy is very likely to 
occur, which is the Jahn-Teller effect. 

The phononic degrees of freedom are very often treated 
as classical variables. The atoms then have well defined po- 
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sitions described by a vector in real space. This assump- 
tion is well justified if the motion of the atoms is very slow 
compared to the spin dynamics, or more precisely, when 
the zero point energy is small compared to all the other 
energy scales. Going back to the spin-Peierls instability in 
the antiferromagnetic S = | spins chain, this assumption 
will fail if the change in energy due to the distortion is 
small compared to the zero point energy of the phonons. 
In other words, if the dimerization is small, the zero point 
motion will restore an undistorted ground state as was 
shown by a quantum treatment of the phonons |5l(il7l8| . 

In the following, we are interested in small magnetic 
rings. This kind of systems in the presence of spin-phonon 
coupling was already studied with a classical treatment 
of the phonons |3] In this paper, we pay special atten- 
tion to the quantum nature of the phonons and show that 
some important departures from the classical picture may 
occur, depending on the different parameters which are 
the phonons frequency, the exchange interactions and the 
spin-phonon coupling. We will concentrate on an antiferro- 
magnetic 4-site cluster and study it in great details, treat- 
ing the phonons in a fully quantum mechanical way by 
means of exact diagonalizations of the hamiltonian. We 
also restrict ourselves to a "weak" spin-phonon coupling 
since this is the most relevant experimentally, but consider 
all possible values of the two energy scales involved in the 
problem : the exchange interactions (J) and the phonons 
frequency (uj). We also present a classical treatment of 
the phonons and show that this approximation is valid in 



2 



Maged Elhajal, Frederic Mila: Dynamic spin Jahn- Teller effect in small magnetic clusters 



the limit J> hw. In this limit, a spontaneous distortion 
is favoured since the gain in magnetic energy (oc J) is 
much larger than the cost in elastic energy (oc Hlo) and 
the ground state is two-fold degenerate. Taking into ac- 
count the quantum nature of the phonons will cause some 
tunneling between these two classically degenerate ground 
states and lift the degeneracy. This is in complete analogy 
with the dynamic Jahn- Teller effect in electronic systems 
U, although in our case the two-fold degeneracy is not 
due to a degeneracy of the spin system but results from 
the static distortion. A semi-classical approach taking into 
account the tunneling between the two distorted states 
will be presented and the results will be compared with 
exact diagonalization to show quantitatively the limit of 
this method. In the opposite (anti-adiabatic) limit, the 
phonons energy is much higher than the exchange inter- 
action and this is then the domain of validity of the pure 
Heisenberg model at low energy. If the phonons energy is 
lowered, but still higher than the magnetic interations, a 
scheme to integrate out the phonons degrees of freedom |U] 
is very often applied, which results in additionnal effective 
exchange interactions (renormalization of the exchange 
constants and appearance of effective longer range inter- 
actions) . It is shown that its usefulness for small magnetic 
clusters is limited. The (smooth) transition from the adi- 
abatic (J 3> fno) to the anti-adiabatic (J <C fajS) regimes 
is shown on the basis of exact diagonalization. Finally, 
a possible interpretation of the low energy spectrum of 
Cu 2 Te 2 05Br2 is presented. 



2 Model 

Starting from a Heisenberg model: 



the coupling to the phononic degrees of freedom is intro- 
duced as a dependence of the coupling constants Jy (dy ) 
on the distances dij between the sites i and j. This is al- 
ready a simplification of what a realistic coupling may be 
since we assume that Jjy depends only on the distance be- 
tween the atoms and not on their exact relative positions 
(Rj — Rj) as it does in reality. We are however interested 
in general trends of what happens in the presence of spin- 
phonon coupling and not in a fully realistic model. For the 
same reason we assume the simplest model for the motion 
of the atoms : They form independent harmonic oscillators 
of frequency lu. This is Einstein's model of dispersionless 
phonons. 

In transition metal ions, the magnetism is due to d 
orbitals and the magnetic interactions are due to overlaps 
of these orbitals. In Anderson super-exchange model [TH] 
(relevant for Mott insulators) the exchange interaction is 
J oc jj, where t is the hopping between magnetic sites 
and U the on-site Coulomb energy. The dependence of 
J on the distance d between the atoms comes through 
the dependence of the hopping t on d. It can be shown 
[TT] that for d orbitals t ~ i which leads to an explicit 



expression for the spin-phonon coupling : J(d) = Jo (^j) , 
where do is the distance between the atoms at rest and a a 
number whose value is 10 according to the previous simple 
arguments, but whose actual value in real compounds may 
be slightly different. If d deviates only slightly from do an 
expansion of the previous expression leads to : J(d) = 

Jo (l — a d ~£ a ^ which will be used in the following. 

The behaviour of a cluster of spins in the presence of 
spin-phonon coupling depends on the number of sites and 
on the spin value [3]. However, this paper concentrates 
on the effects of the quantum nature of the phonons and 
for this purpose, we restrict ourselves to a ring of 4 spins 
where the atoms can only move along the ring. Let Xi be 
the deviation of the atom i with respect to its position at 
rest, then the hamiltonian is written as : 



2m 



(i) 



1 - — (x i+ i 
do 



where the first two terms are the kinetic and potential en- 
ergy of the phonons. If the phonons are treated as quan- 
tum oscillators then it is more convenient to make the 
usual transformations to the operators a and a + and ijljl 
becomes : 



H = fau (al a * + \ 

[i - p (oi+i + 4+i 



(2) 



Si+i-Sj 



with 



d 



Ivtilu 



(3) 



In the hamiltonian J2J), a, do and m are microscopic pa- 
rameters and in the following they will be fixed to some 
typical value. Although there are only two energy scales 
(J and hu>), it is not possible to reduce the parameters 
in to only one relevant parameter (for example ^j-) 
because (3 oc ^= as is seen from ©. However, to avoid 
dealing with two independent parameters, in the follow- 
ing lu is fixed to a typical value of = 100-fT which is a 
realistic value for phonons in real compounds, and J will 
be varied. For this value of ui, taking typical values a ~ 10, 
do ~ 3.5 A and m ~60 atomic masses, we find f3 ~ 0.2. 
Except when specified, [3 will be fixed to this typical value 
in the following. For other values of j3 the physical pictures 
presented remain valid although the results may depend 
quantitatively (but not qualitatively) on (3. An example of 
this dependence is investigated in more detail in sectional 
where (3 is varied. 

Strictly speaking, one cannot express all the results as 
a function of the ratio 4-, because (3 depends also on u> and 
the calculations are actually done for one fixed value of lu. 
However, keeping in mind that we have fixed = WOK, 

we will express all the results as a function of ^ because 
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the relative values of these two energy scales is the relevant 
parameter which determines the characteristic behaviour 
of the system. 



3 Exact diagonalizations. 

In order to keep the quantum nature of the phonons we 
will diagonalize exactly the hamiltonian J5J. We first give 
some technical details about exact diagonalizations, which 
can be skipped by the uninterested reader. 



3.1 Technical details 

The dimension of the Hilbert space of a quantum harmonic 
oscillator is infinite and it should be truncated in order to 
be diagonalized numerically. This is done by introducing a 
cut-off in the maximal number of phonons on each site. For 
each set of parameters, we estimated the error introduced 
by this cut-off by inspection of the evolution of a subset of 
energy levels as the cut-off is varied (see figures ^ and [2J . 
For all the results presented in the following sections, the 
error bars would be completely negligible and for this rea- 
son they were omitted. The phonon wave function was de- 
scribed using the basis in which the number of phonons is 
diagonal. This is a good choice in the limit oiftio ^ J since 
the number of phonons is then "almost" a good quantum 
number (see section 01 . The convergence, with respect to 
the number of phonons states taken into account, is then 
very fast because the eigenstates will be correctly approx- 
imated using a small number of vectors from the choosen 
basis. Indeed, for the largest values of fuv considered in 
this work, it was found that taking into account a maxi- 
mum number of 5 phonons per site is already an excellent 
approximation. As hu> is decreased, more and more excited 
levels of the harmonic oscillators are involved. In principle, 
to build the (classical) static distortion found in section^] 
all the energy levels of the harmonic oscillators are needed. 
For the largest values of J considered, up to 40 phonons 
per site were taken into account in order to get a good 
convergence of the results. This represents a huge Hilbert 
space and it was necessary to take into account both geo- 
metrical symmetries of the cluster and the fact that S z is a 
good quantum number in order to be able to diagonalize 
numerically the hamiltonian, using Lanczos's algorithm. 
The diagonalization was performed separately in each of 
the S z — 0, 1 and 2 subspaces (S z = — 1 and — 2 being 
respectively degenerate with S z — 1 and 2). The geomet- 
rical symmetries we took into account are the 4 rotations 
(including identity) around the axis passing through the 
middle of the ring and perpendicular to it. These rota- 
tions perform circular permutations of the 4 sites and if 
the ring is seen as a linear chain with periodic boundary 
conditions then they are equivalent to 4 translations along 
the chain. There are 4 irreducible representations for this 
group of symmetries which we label with k — 0, 1, 2 and 3. 
The diagonalizations are performed separately in each of 
the subspaces associated with these 4 irreducible repre- 
sentations. 
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Fig. 1. Evolution of the low energy spectrum with the maxi- 
mum number of phonons allowed on each site. Symmetries are 
used, and only the lowest energy levels corresponding to the 
irreducible representation (S z = 0, k = 0) are shown. 
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Fig. 2. Detailed view of the third group of levels from figure 
Although the rough convergence seems to be achieved for 
30 phonons per site (and indeed the ground state energy is 
already converged to ~ 10 -6 ), more phonons are needed to get 
the fine structure of the spectrum (notice the fake degeneracy 
lifting for 30 phonons per site). 



3.2 Results 

The evolution of the spectrum as a function of J- obtained 
from exact diagonalizations is shown in figure 1 ^] Several 
aspects of this spectrum will be studied in more details 
in the next sections and only a few remarks will be made 
here. 

First of all, for small -i- where the phonons are es- 
sentially frozen at low energy, one can recognize the four 
levels of the pure Heisenberg model (without coupling to 
the phonons). In the opposit limit (large it is seen that 
the ground state is two-fold degenerate, which is the quan- 
tum equivalent of the static distortion predicted when the 
phonons are treated classically (see figureQJ- However, the 
value of -i- at which the distortion appears in the classical 
approximation (arrow on figure [3J| is quite different from 
that at which the degeneracy of the ground state appears 
(see section 0] for more details). The continuous evolution 
of the spectrum from one limit to the other is one of the 
original results of our work. Several methods exist to solve 
the problem starting from one limit or the other. They will 
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j/huj 

Fig. 3. Transition from the anti-adiabatic (hu S> J) to the 
adiabatic (hu <C J) regimes. For each value of J- only the 
lowest states are represented. In the anti-adiabatic limit, the 
low-energy spectrum is made of the four levels of the pure 
Heisenberg hamiltonian. As J- increases, initially higher en- 
ergy levels involving phonon degrees of freedom come into play. 
For large values of the system enters the adiabatic regime 
where the phonons behave classically. The ground state is then 
two-fold degenerate, since one singlet state (empty squares) 
falls on the ground state after an anti-crossing with a level of 
higher energy (full squares). For comparison, the critical value 
of jr- above which a static dimerization is predicted in a clas- 
sical treatment of the phonons is shown. 



be presented in the next sections, which will help under- 
stand the spectrum of figure|21and the domains of validity 
of these approximate methods will be made explicit by 
comparison with the results of exact diagonalizations. 



4 Classical (adiabatic) and semi-classical 
phonons 

Small rings of spins with adiabatic phonons were studied 
and several interesting features were found by varying the 
number of sites and the spin . The calculation for the 4- 
site S — | ring is repeated here for the sake of comparison 
between the classical (adiabatic), semi-classical and fully 
quantum approaches in the next paragraphs. 



4.1 Classical phonons 

The ground state of the hamiltonian (JTJ is achieved when 
the atoms are at rest (pj = 0) and at the minimum of the 
total energy : 



E : 



IE*? 



where Em is the energy of the spin system ground state for 
fixed values of the positions Xj. Em is calculated quantum- 
mechanically. Minimizing this expression with respect to 
the Xj yields the energy of the ground state : 



E 



-J- JW1 + 3 



2a V 
*7 



do + 2x 




do — 2x 

Fig. 4. Two degenerate ground states obtained with classical 
(adiabatic) phonons if the spin-phonon coupling is above the 
critical value (equation [IJ. 



with x the displacement as represented on figure This 
expression has a minimum at finite x only for large enough 
spin-phonon coupling constants. The critical value is : 



a c = d o\I^J 



(4) 




For a < a c the atoms will remain undisplaced, whereas 
for a > a c a static distortion, as represented on figure 0] 
will occur. The amplitude of the distortion is then : 



(5) 



4.2 Semi-classical approach 

The preceeding classical method neglects the quantum 
fluctuations of the phonons and it can be improved by 
reincorporating them on top of the classical solution. In 
order to do this, we first solve the quantum problem in a 
mean-field approximation inspired by the classical dimer- 
ization of figure 0] This is done by introducing "by hand" 
a static distortion : 

8 = {ai+af) = (a 3 + a^) = -(a 2 + a^) = -(04 + 04) (6) 

which plays the role of the classical variable x in the 
classical treatment since for a harmonic oscillator (x) = 



2kx 2 



77~ — (a + a + ). The two distorted states of figure Ulcor- 

respond to 8 > and 8 < 0. Using JBJ, we replace in the 
hamiltonian the and af by their expression in term 
of 8. Doing so, one obtains a hamiltonian with only spin 
variables which can be diagonalized exactly, and the mean 
values of (Sj.Sj_|_i) in the ground state is calculated as a 
function of 8. Replacing the S^.S^+i by their mean value 
in the hamiltonian J2J), one obtains then a hamiltonian 
which depends on 8 with only phonon variables. This is 
a displaced harmonic oscillator hamiltonian which is eas- 
ily diagonalized and finally, self-consistency is obtained 
by requiring that this phononic hamiltonian leads to the 
same distortion 8 as the one first introduced in the spin 
hamiltonian. 

Not surprisingly, the results arc similar to the adiabatic 
phonons results. The distortion <5 and the critical value (3 C ■ 
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<y = ± 




(3 C 



(7) 



are the same as ijtjj and (0J given the relationship between 
a and (3 (see equation Q) and between 5 and x (see dis- 
cussion after equation (6)). However, unlike in the pure 
classical approach, we have now an explicit wave function 
of the ground state which will be a direct product of a 
phonon part and a spin wave function. The phonon part 
is the wave function of a displaced harmonic oscillator, 
with the displacement given by equation Q. The spin 
wave function is : 



i= v / TTT(i rm) + imt))+ 



Y) (| TTU) 



where Y — 



-^l+^^VT^TKl TUT) 



2V2 



UTT»- 

uu» 



g%p and the two possibilities (±) correspond 



to positive and negative <5 (the two solutions of figure 0}. 

This mean-field solution has neglected the quantum 
fluctuations, and one can improve it by re- introducing 
them, on top of the two-fold degenerate distorted mean- 
field ground states. Technically, this is achieved by re- 
stricting the Hilbert space to a two-dimensionnal space 
generated by the two solutions obtained above. Rewriting 
the full hamiltonian J5J in this subspace, and diagonaliz- 
ing it, one obtains two energy levels : 



E = 
Ei = 



a - 



1 + (1|2> 
a — b 



1 



(1|2> 



where 



J / 1 
a = 2tw- -\2 + Y + — 

(l|i/|2^-^(3 + 2F-r 2 ) 
(1|2) = Ye- 26 ' 



TYhw 



(8) 
(9) 

(10) 

- 2g2 (11) 
(12) 



(1|2) being the overlap between the two mean- field so- 
lutions (the exponential comes from the overlap of the 
phonon wave functions and Y from the overlap of the spin 
wave functions) . It is useful for the following discussion to 
note that b is the matrix element of the hamiltonian con- 
necting the two states and the splitting between them is : 



2Je- 2S ' 2 (l 



Y 2 ) 



l-Y 2 t 



-452 



(13) 



4.3 Discussion 



It is easily seen on the expression {7J) of the distortion 5 
that for large values of the ratio the distortion will be 
important (which is intuitive since the gain in magnetic 
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Fig. 5. Semi-classical approach from the adiabatic limit, (top) 
Comparison of the energies of the two lowest levels (equations 
||SJ| and @) as obtained from the semi-classical approach (solid 
line) and exact diagonalization (points). For larger values of J- 
the energies obtained with both methods will become identi- 
cal (not shown here for clarity). The third energy level from 
exact diagonalization was included, (bottom) Comparison of 
the splitting between the two lowest energy levels from the 
semi-classical approach and exact diagonalization. Although 
the agreement is not excellent, the same exponential behaviour 
appears for large J- in both approaches. 



energy is large while the elastic cost is small). As ■J- and 
S increase, the overlap Ijl2(l and the matrix element (|1 If) 
between the two states will vanish exponentially. The two 
mean-field solutions (identical to the classical result) will 
thus become assymptotically eigenstates of the hamilto- 



oo. This is how the classical behaviour of 
the phonons is recovered starting from the quantum hamil- 
tonian. As is lowered, the matrix element between the 
two states will become more and more important, giving 
rise to a tunnelling probability between the two classical 
states and a splitting will appear. 

The two energy levels obtained, as well as the split- 
ting between the two states, are represented on figure 03 
and compared with exact diagonalizations. As seen from 
the figure, the approach is qualitatively correct since for 
large -^j, the exact spectrum has two energy levels well 
separated from the higher levels and which show the same 
behaviour as predicted by the semi-classical approach. The 
semi-classical approach however fails to predict quantita- 
tively the value of the splitting although it has the correct 
exponential behaviour at large 4-, 

Finally, in the regime where J- is very large, a very 
simple physical picture arises. The system is expected to 
be classically dimerized and the low energy degrees of free- 
dom are the vibrations around the distorted states. They 
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Fig. 6. Low energy spectrum obtained with exact diagonal- 
ization in the adiabatic regime (J ^> hui). All states are singlet 
states, triplet states being much higher in energy. The ground 
state is two fold degenerate, which is the quantum equivalent 
of the two classically dimerized states. The rough structure of 
the excited states correspond to harmonic vibrations around 
the dimerized states. The splitting within one group of close 
lying states can be seen as a perturbation of these oscillations 
by the magnetic degrees of freedom. 



are thus displaced harmonic oscillators with the same fre- 
quency as the phonons. This scenario is confirmed by exact 
diagonalizations (see figure 



5 Expansion from the anti-adiabatic limit 

In the abscence of spin-phonon coupling the eigenstates 
of the hamiltonian are direct products of a magnetic and 
phononic parts. Furthermore, if fttu 3> J, the energy levels 
will be grouped around the energy levels of the harmonic 
oscillators (which are the large energy scale) and within 
each group there will be a fine structure given by the mag- 
netic part. The lowest group of states is nothing else but 
the spectrum of the pure Heisenberg model, since these 
states have no phonon. Including spin-phonon coupling, 
one can consider this term as a perturbation compared to 
the pure Heisenberg and the pure phononic terms. This is 
because (3 is a rather small number and because the low- 
est (unperturbed) states have zero phonons and thus the 
phononic operators at and af will give contributions of 
the order of to second order in perturbation. Restrict- 
ing the attention to these zero phonon states, one can take 
into account the spin-phonon coupling in perturbation and 
interpret the result as an effective spin hamiltonian, thus 
integrating out the phonons degrees of freedom [£f[. To 
second order in perturbation theory, this will introduce 
effective second neighbours interactions and renormalize 



E/J 



-B B B B— 



0.2 0.4 0.6 0.8 1 

J/TlUJ 

Fig. 7. Low energy spectrum as a function of -S- . The lines are 
the spectrum obtained with a pure Heisenberg model renor- 
malized by high energy phonons to second order in pertur- 
bation (see text and equation 1141 '). The points are obtained 
from exact diagonalization of the hamiltonian (|5J. For clarity 
only the lowest energy levels were represented at each value of 
J/huu, and an energy of 2hu> (zero point energy of the 4 phonon 
modes) was substracted. S=0 states are represented by circles, 
S=l states by triangles and S=2 states by squares. The anti- 
crossing between two S=0 states (full circles) is emphasized. 
The lowest of these two states will be the first excitation be- 
low the first triplet state for ^ > 0.9. In the adiabatic limit 
(J S> huj) it will fall onto the ground state (see figure|2J. 



the first neighbours interactions, so that the total hamil- 
tonian has the usual Heisenberg part and in addition : 



(14) 



Two distinct assumptions must be fulfilled in order for this 
effective hamiltonian to be valid, although they are both 
related to the fact that Hu> ^S> J. The first one is the usual 
validity of perturbation theory, and one could improve the 
above scheme by going to higher orders in perturbation. 
The second assumption is that the phononic degrees of 
freedom (which are integrated out) remain at higher en- 
ergies than the spectrum obtained with the above effec- 
tive hamiltonian. In other words, the states which have 
one or more phonons should remain well separated from 
the lowest energy levels for which we derived the effective 
hamiltonian. 

Exact diagonalization results are compared with the 
energy levels of the pure Heisenberg hamiltonian, renor- 
malized by l|14|) (see figure 0). For very small the low 
energy spectrum is essentially that of a pure Heisenberg 
hamiltonian (4 initial levels for small and the per- 
turbative approach reproduces well the small departures 
(~ 1%) as 4- increases. For ^ ~ 0.4, states from higher 
energy involving phonons will fall on the lowest levels, 
causing great departures from the (renormalized) Heisen- 
berg model. The perturbative approach fails completely 
to reproduce such a change in the spectrum since it inte- 
grates out the phonons. This commonly used method (for 
other systems) appears thus useless for small magnetic 
systems since, in its range of validity it only predicts de- 
partures of ~ 1% in the renormalized exchange constants, 
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which is orders of magnitude smaller than the perturba- 
tion brought by higher energy phonons as seen in exact 
diagonalizations. 



6 Physical applications 

Appart from testing the validity of various approximate 
methods, two aspects should be underlined from the re- 
sults of the previous sections, which are important for 
small magnetic clusters, such as molecular magnets. First, 
the spin-phonon coupling can tremendously change the 
low energy spectrum of small magnetic clusters. This was 
already pointed out by a previous work where the phonons 
were treated classically Second, the quantum fluctua- 
tions associated to the phononic degrees of freedom must 
be taken into account to describe correctly the low energy 
spectrum, except for the classical limit (J ^> fkS), which 
is unlikely to occur in molecular magnets, given the small- 
ness of the exchange interactions in these systems, which 
are typically much smaller than the phonon frequencies. 
It was also shown that the spin-phonon degrees of free- 
dom will qualitatively change the low energy spectrum in 
small magnetic clusters, even in the abscence of a static 
distortion, a phenomenon which is similar to the dynamic 
Jahn- Teller effect for the electronic degrees of freedom. 

For instance, it was shown how a spin-phonon cou- 
pling can be responsible for a singlet-singlet first excita- 
tion in a 4-site cluster (whereas the lowest transition is 
a singlet-triplet one in the absence of spin-phonon cou- 
pling). The value of the spin-phonon coupling was so far 
fixed to p = 0.2 and to obtain a singlet first excitation 
requires jtj > 0.9. Higher values of the spin-phonon cou- 
pling would however require lower values of -d- which is 
more likely to occur in real materials. In order to make 
this argument more quantitative, figure [3] shows the de- 
pendence on J- of Pss, the critical value of P needed for 
the first excited singlet to be lower in energy than the first 
triplet. As can be seen from this figure, even if the phonon 
frequencies are up to two times the exchange interactions, 
reasonable values of P are expected to be enough to change 
qualitatively the spectrum and to bring a singlet excita- 
tion below the first triplet. 

The magnetic properties of the compound Cu2Te205Br2 
are not fully understood. Although magnetic measure- 
ments ^2] and ab-initio calculations 13 point to a system 
of weakly coupled plaquettes of 4 spin--| copper ions, a 
magnetic order was recently observed |lo| . Partial expla- 
nations were given |14I15| but an overall understanding of 
all the features is still lacking, the main puzzle being the 
existence of both magnetic order and a low lying singlet 
excitation. A coupling between the plaquettes is needed 
for the appearance of magnetic order and the observed 
phase transition. However, some features, such as the ex- 
istence of a low lying singlet excitation can be understood 
on the level of one isolated plaquette. 

Experimentally, Raman spectroscopy showed that 
the ground state as well as the first excitation are sin- 
glet states. Due to the symmetry of the cristal, only two 




Jl 
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Fig. 8. Elementary plaquette of four spin-i Cu 2+ ions in 
Cu2Te20sBr2 with first (Ji) and second J2 neighbours cou- 
plings, ab-initio calculations find negligible J2, whereas Raman 
spectroscopy shows a first singlet excitation which implies sig- 
nificant J2 in a pure Heisenberg model. It is proposed that 
indeed J2 is negligible and the first singlet excitation is due to 
spin-phonon coupling 
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Fig. 9. Critical value of the spin-phonon coupling /3 needed 
in order to bring the first singlet excited state below the first 
triplet. 



coupling constants can appear in one plaquette, as repre- 
sented on figure 03 If Ji is of the order of Ji then the 
ground state, as well as the first excited state will be 
singlet (S=0) states. However, ab-initio calculations have 
shown that J2 is in fact negligible ^Sj) in which case the 
first excitation is a triplet state and the first singlet ex- 
citation is twice higher in energy, in contradiction with 
experimental results. We propose that spin-phonon cou- 
pling may be responsible for a first singlet excitation, in 
the abscence of second nearest neighbours coupling (J2). 
Figure |3 indeed shows that for J > hu, the first excited 
states are two almost degenerate singlet states, below the 
first triplet state. The singlet excitation is actually not 
a pure magnetic excitation, as would be obtained by a 
next-nearest neighbour coupling J2, since it involves also 
phononic degrees of freedom. Although the present model 
is over-simplified for a quantitative description of the com- 
pound Cu2Te205Br2, it shows that such a mechanism is 
possible. 



7 Conclusions 

We have studied the effect of spin-phonon coupling in a 
4-site S = \ antiferromagnetic ring taking into account 
exactly the quantum nature of the phonons by means of 
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exact diagonalizations. Perturbative and variationnal ap- 
proaches in the adiabatic and anti-adiabatic limits were 
presented and compared to exact diagonalizations and a 
quantitative estimate of their domain of validity and ac- 
curacy was presented. The evolution of the spectrum from 
the classical regime to the anti-adiabatic one was shown 
by means of exact diagonalizations. From these results, it 
is expected that spin-phonon coupling might play an im- 
portant role in small magnetic clusters, as for instance in 
molecular magnets. Furthermore, it is very likely that tak- 
ing into account the quantum fluctuations of the phonons 
is necessary to achieve a correct description of these sys- 
tems. This is examplified in the magnetic analog of the 
dynamic Jahn- Teller effect which was underlined. An in- 
terpretation of the low energy spectrum of Cu2Te205Br2 
was proposed. 

We have concentrated on a 4-site ring, but we believe 
the results are valid for other clusters, at least qualita- 
tively. This is supported by the similarity of even clusters 
if the phonons are treated classically 3 . There might how- 
ever exist a discrepancy between clusters with even and 
odd numbers of sites. Indeed, even-rings have a 2 fold de- 
generate distorted ground state, whereas this pattern is 
frustrated in the case of odd clusters, leading to a much 
higher (classical) degeneracy. The quantum fluctuations 
might thus play an even more important role in odd clus- 
ters. Another significant parameter could be the value of 
the spin, since this dependence was observed in clusters 
with a classical treatment of the phonons • These prob- 
lems are left for future investigation. 
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